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Abstract 

Peierls distortion and quantum solitons are two hallmarks of 1-dimensional condensed- 
matter systems. Here we propose a quantum model for a one-dimensional system of non- 
linear ly interacting electrons and phonons, where the phonons are represented via coherent 
states. This model permits a unified description of Peierls distortion and quantum solitons. 
The non-linear electron-phonon interaction and the resulting deformed symmetry of the 
Hamiltonian are distinctive features of the model, of which that of Su, Schrieffer and 
Heeger can be regarded as a special case. 

One-dimensional condensed-matter systems have attracted increasing inter- 
est in several branches of physics: not only do they have promising appli- 
cations in information-processing technologies, but they also play a central 
role in biological molecules. They consist of linear chains of 'ions' whose 
conduction electrons move primarily along the chain axis. Hence, attention 
can be confined to a single chain. 

In this paper we shall consider a half-filled 1-d chain, i.e., one consisting of 
an even number of ions, each of which carries a single conduction electron. 
This class of systems is particularly interesting because they undergo Peierls 
distortion: at equilibrium the ions shift from the equally-spaced configura- 
tion and assume a dymerised pattern, where the bonds between adjacent 
ions are alternatively short and long. To date, the best account of this phe- 
nomenon is Peierls's theorem, [T], stating that the dimerised configuration 
minimises the total energy of a 1-d half-filled chain. However, the proofs 
of the theorem rely on models which represent the ion coordinates as static 
variational parameters, |2j,[3 |. Hence, they do not clarify how Peierls dis- 
tortion affects the energy spectrum of the whole system, let alone how it 
contributes to collective, dynamical effects that are also expected to arise. 
Specifically, since the equally-spaced configuration of identical ions has a re- 
flection symmetry, there are two topologically non-equivalent ground states 




Figure 1: Peierls distortion. The unstable configuration, E; the two stable vauca, 
V\ and V2; a soliton S interpolating between them. 

(vacua) in which Peierls distortion may result, one obtainable from the other 
by exchanging the position of long and short bonds (fig [I]). Hence, there are 
additional stable states of the system, known as quantum solitons, where 
the two degenerate vacua coexist (at a given time) and a kink-like domain 
wall (S in fig{T]) interpolates between them. Soliton peculiar properties, such 
as fractional charge eigenvalues, have thus far been described only by phe- 
nomenological models [1], [S], [BJ, and so has their connection with Peierls 
distortion. In fact, the best available model addressing the latter issue — 
proposed by Su Schrieffer and Heeger (SSH) for polyacetylene [7]— assumes, 
for phenomenological reasons, the electron-phonon coupling to be linear in 
the ion displacements, as the expected values of displacements in the dis- 
torted ground state are much smaller than the interatomic distance. Overall, 
we lack of a unified description for the systems sustaining Peierls distortion: 
the available one is disconnected and incomplete, since it relies on ad hoc 
models, valid only in special regimes. 

In this paper we propose a second-quantized model for a 1-d system support- 
ing Peierls distortion, where the electrons and phonons interact non-linearly. 
The ion-coordinates are described as semiclassical dynamical variables in co- 
herent states. In this framework we prove Peierls's theorem, showing that 
Peierls distortion is a direct manifestation of the non-linear electron-phonon 
interaction, and we show that the system supports a kink-like excitation 
propagating along the chain at constant energy. 

The model Hamiltonian and its quantum symmetry — Our model describes 

a linear chain of 2L ions interacting with 2L conduction electrons. 

2L-1 

The Hamiltonian is: H = H p h + H e i_ p h. Here, H p h = | (p| + itj) is the 

j=0 
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Einstein-phonon Hamiltonian with the ion mass, oscillator frequency and h 
set equal to 1; Uj is the ion displacement from the equilibrium position Rj 
(see fig. [TJ and pj is the corresponding momentum, {\uj,pk\ = ihSj^)- The 
degrees of freedom in the plane orthogonal to the chain are considered as 
'frozen'. The electron-phonon Hamiltonian, 

2L-1 

Hei~ph = ~ ^2 • 1 Uj ■ : li + h - c ) 

3=0 

is a tight-binding Hubbard Hamiltonian written in terms of fermionic cre- 
ation and annihilation operators /j, fj ({fj,fj} = o~j.e, {fj,fe} = 0). The 
index j includes both electron position and spin, but the latter is irrelevant 
in the present context and will not be explicitly written. The 'amplitude ' 
for an electron 'hopping' from site j to site j + 1 is the operator tjj+i, 
depending on the phonon degrees of freedom. It includes, as argued in [5], 
the electron-phonon coupling: 

hj+l = T ex P (C (uj+i - uj)) exp (k (pj+i - pj)) C, k G M, T G C, 

where £ and k depend on the form of the Wannier functions out of the ion 
core [9]. Even though no explicit phonon-phonon interaction is considered, 
the phonons are indirectly coupled via H e i_ p h. The SSH Hamiltonian can 
be recovered by setting k = and approximating tjj+i to first order in £. 
The Hamiltonian H e i_ p h has a local dynamical symmetry and a global 
symmetry, both associated with the Hopf quantum algebra U g (su(2)), [5]. 
The generators of this algebra {K^\ K^} close the g-deformed commu- 
tation relations, [10J, [kW,K&] = ±K^\ [K(+\K^] = [K^] q , where 

A — A 

[A] q = 9 q Zg-i an d Q is the deformation parameter (which can be assumed 
to be real). The algebra U q {su{2)) belongs to the universal envelope of su(2) 
and reduces to the latter for q — > 1. As discussed in [8], the generators of 
the local dynamical symmetry contain both fermionic and bosonic operators, 
and so do the generators of the global symmetry (defined via co-producting 
the local operators). Also, q depends on the physical constants of the model, 
in such a way that these symmetries reduce to a non-deformed su(2) symme- 
try (i.e., q — > 1) if £ — Y or k — >■ 0. Since both ( and k are non-negligible 
in the systems supporting Peierls distortion, we expect that the quantum 
symmetry will play a central role in describing this phenomenon. 
Staggering in phonon coherent states — In view of the different time-scales 
of the phonon and the electron dynamics we adopt for the phonon operators 
a semiclassical, dynamical description in terms of Glauber coherent states 
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[TT] . The global coherent state of the ions is \Z) = (^) \zj), 1 < j < 2L — 1, 

i 

where \zj) = e~ 5 I z jI e 2ja ? \Q)j, dj is the j-th single-phonon creation operator, 
1 0) j is the vacuum state such that Oj|0)j = and Zj £ C. Here, Rezj and 
Imzj represent respectively the expected values in coherent states of the j-th 
ion displacement and momentum. 

Inspired by [7], we first explore the possibility of a staggered ground state 
by setting Zj = {—) 3 z, where z is a variational parameter to be found min- 
imising the ground-state energy (z = corresponds to the non-dimerised 
configuration, E in fig|T]). The staggering condition is well defined, as it 
involves the expected values of the ion positions and momenta in coherent 
states; in constrast, in [7J the condition is imposed on the ion position op- 
erators, ignoring the effect on the conjugate momenta, which may lead to 
a difficulty with Heisenberg principle. Besides, global momentum conserva- 
tion implies, in our semiclassical picture, that the ion momenta are staggered 
too. We also require T = texp(—ih(^K), with t S K, so that (Z\tjj+i\Z) G R 
(time-reversal symmetry) . 

The averaged Hamiltonian and its symmetry — Representing the Hamilto- 
nian in coherent states we find: H = (Z\H\Z) = H p h + Hei-ph = {%\H p h\Z) 
+(Z\H e i_ ph \Z) where H ph = 2L (4(Rez) 2 + (Imz) 2 + |) (a c-number) and 

2L-1 

Hei-ph = ~9^2 ( cosh (3) ~ (~y sinh (3)) (fj+ifj + h - c -) . 

j=0 

where 'state location' 3 = 2y/2 (Re(C,z) + Im (kz)) and effective coupling 
g = texp (£ 2 + k 2 ) have been introduced. 

Fourier-transforming the {fj} into the standard particle- hole fermionic op- 

L-l 

erators {c^, v^}, < k < L — 1, we have H e i- P h = Hki with 

fc=o 

H k = -2eJ 3 - 5 [J+ + J_] , (1) 
where e = e(k, z) = g cosh(3) cos (j;k), 5 = 5(k, z) = g s'mh.(^) sin (j;k) and 
J + = v\ck, J- = j\, Jz = \ (j 1 ^ ~ n ^^j (dropping the mode-index k for 
simplicity) . 

The operators J a ,a£ {+, — , 3}, close an su(2) algebra in the spin-^ repre- 
sentation, T>i. H e i_„h has therefore the dynamical symmetry described by 

2 y 

the algebra A = (Bfc su (2)(fc)> whereas the original Hamiltonian had both 
a global symmetry and a local dynamical symmetry associated with the 
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q— deformed algebra U q (su(2)). Indeed, the generators of the quantum sym- 
metry [S], when averaged in coherent states, lose their dependence on the 
phonon operators and reduce to the generators of su(2). Since the quantum 
symmetry is induced by the electron-phonon interaction and the latter is 
central to the description of systems supporting Peierls distortion, we shall 
now define a procedure to restore it. 

Restoring the quantum symmetry — To this end, it is not convenient to refer 
to Uq(su(2)), as it is not a proper algebra and does not provide the group 
operation to diagonalise the Hamiltonian. Instead, we define a proper al- 
gebraic structure, A q , having more appealing properties. Specifically, A q 
is the three-dimensional submodule of U q (su(2)) closed with respect to the 
deformed adjoint action, [ , ] : [K^\f] g = K^fq K(3) - q^^fK^, 

[K^ ,g] q = K®g - g K® and [fg, h] g = [f, [g,h ) q ) q , V/, g,heU q («*(2)) . 
Resorting to the deformed adjoint action is necessary since there is no three 
dimensional submodule closed with respect to the deformed commutation 
relations in U q {su{2)). A q is generated by the operators 

where are the generators of U q (su(2)), £ q = 2q w (q + q x ) 

and g,ui£i Also, = . Notice that for q — >• 1 A q coincides with 

su(2), as [ , ] x = [ , ]. 

In order to restore the quantum symmetry we simply replace in each 
operator J a with Ja \ for every k, i.e., with 

= -2e4 q) - 6 + J W ) . 

The numbers q and w therefore become parameters of the model. This 
Hamiltonian, formally identical with ([!]), is endowed with a quantum dy- 
namical symmetry which retains memory of the of the original Hamilto- 
nian symmetry. Indeed, as a side remark, a physical interpretation of the 
quantum symmetry restoration may be provided by writing the generators 
of A q in terms of "g-deformed" fermionic-like operators, f, ft, obeying de- 
formed Clifford anticommutation relations: {f , ft} = — f 2 = ft 2 = 

— \ Q ^ sm g ^ ^ , Q ~ log(g). These operators represent "dressed " electrons, 

retaining some memory of the interaction with phonons. 

In the Di representation, A q ~ u(2) = u(l) © su(2) and the eigenvalues of 
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are: 



A± 



e(q - q- 1 )^ =F Jq 2 "e> + H q 5 2 = A±(fc, , 



(which, in the limit q — > 1, tend to the eigenvalues of H^ ) Since in Di the 
co-product is primitive (A q ~ u(2)) the global symmetry is automatically 
restored by setting H e i_ p f l = . The energy spectrum of the system 

k 

is therefore the sum over the modes k of the eigenvalues A±. 

Proving Peierls's theorem — The ground state energy density has the form 

1 L-l 

£ q {z) = — A+(fc, z) (taking into account the spin degeneracy by a factor 



L 

k=0 

2). In the limit L > 1, 



{*) / 

J 



&k\ l — m q sin (k) 



with p q (z) = - gq w cosh(3), m q = (1 - y 2 ), y q = y 7 ^ tanh(3). 
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Figure 2: The ground state energy density when 1 < £ 9 < 2 and q = 1.5. 

The integral in £ q {z) for < £ q < 1, i.e., < m q < 1 is an elliptic integral 
of the second kind; whereas for 1 < £ 9 < 2, — 1 < m g < 1, and in general 
for £ q > 2, it is the hypergeometric function 2-^1(5, — \, 1; m g ). The integral 
is real and converges if |j| < Cm, where Cm is the value for which \m q \ = 1 



[T2] . This allows one to perform a detailed study of the total energy density 
of the ground state, Eg.s. = £ph{z) + £ q {z). 

For appropriate values of q and £ g (i.e., w) Sg.s. exhibits (see fig. ([2])) a saddle 
point in Hez = = Imz (corresponding to the equally-spaced configuration) 
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and two degenerate minima in ±(Rez,Imz) 7^ 0, corresponding to the two 
degenerate ground states induced by Peierls distortion. This shows that 
the dimerized configuration minimizes the total energy and proves Peierls's 
theorem, as promised. For q — > 1 the total energy has only one maximum 
with zero second derivative, describing a marginally stable equilibrium for 
z = 0. This limit includes, for k = and 3 <C 1, the SSH model (linear 
in Re(z)). This shows that the deformed symmetry is decisive to describe 
Peierls effect. 

Dynamical equations in the ground state — In the coherent-state formalism 
x = 2Rez and p = 2Imz are dynamical variables specifying the ion canonical 
coordinates. £g.sX x >p) can then be considered as the (classical) Hamiltonian 
describing the dynamics of the system ground state in phase space C 2 [15] . 
The motion of the representative point (x(t),p(t)) describes the collective 
evolution of the phonons interacting with the electrons in the ground state. 
The corresponding equations of motion along the line x = p is (see fig(2]) 
represent a non-linearly damped, non-linearly driven oscillator: 



(x -V)(l- V x ) - xV x 



where V x = and 



7T Zgiq + q- 1 ) \ " qJ ^m q [cosh(V2(x + b))]y 

with E(m q ) = 2^1(5, —5, l;rn q ) and F{m q ) = 2 -Fi(|, |, l;m 5 ). The pres- 
ence of a damping factor shows that the Peierls-distorted ground state is 
robust against perturbations. Besides, this non-linear equation strongly 
suggests the existence of soliton excitations. 

The existence of solitons — In analogy with the procedure used to prove 
Peierls's theorem, we explore the possibility of an excited kink-like state at 
site n. To do so, we adopt as a trial description of the phonon coherent state 

n N-1 

\2 n ) = (^) \zj) (^) \zi), where the deformed staggering condition 

j=0 i=n+l 

z - = (-y z , for < j < n , z t = {-f +1 z , for n + 1 < i < N - 1 , 
represents the presence of a kink at site n. The averaged Hamiltonian H' n = 

{Zn\H e l—ph\Zn) IS 

N-1 N-1 

H'n=Y, "ifj+lfi + "nfUjn ~ 2 S £ i )'./}.,./) + h.C. , 
j=l j=n+l 
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where oj£ = g — (c + (—) e s), s = gsinh^), c = gcosh($), and the energy 
spectrum can be obtained by repeating the diagonalisation procedure via 
pseudo-fermionic operators. To prove the existence of solitons we shall argue 
that the dynamics generated by H' n , for appropriate initial states, induce 
the kink to move spontaneously from site n to site n + 1. 
The state at time t, consisting of the fermionic (pseudo-spin) and bosonic 
components, can be written as \ip(t)) = e~2l 2 WI 2 e 2 W at \O)0e~ ttH ' n ^ z ^ |s(0)), 
where z(t) satisfies the canonical equations z = §| , z = — §| and £(z,n) 
is the lowest eigenvalue of H' n , representing the first excited state of the 
system. If the time 5t of the kink motion from site n to n + 1 is very small, 

one has \ip{t + 8t)) = 1 + 6t (i(at - z) <g> I - il <g> (H' n + z ^gf)) |^(t)) to 
first order in 5t and ^f- ~ ||P^ +1 ' n \ where 

P (n+l,n) ^ x> n+1 -H' n = UJ n {fl +2 f n+l - f n+1 f n + h.C.) . 

This operator has a doubly-degenerate O-eigenvalue, whose corresponding 
eigenspace, Vo, is spanned by the states ■ ■ ■ ® \n n ) ® \n n+ i) <8> \n n+ 2) • • • , 
with rij = 0, 1, and (conserved) total number n ra + n n+ i + n n+ 2 = 1. 
Suppose now that IV'(i)) ^ s a superposition of {0^. n ' > }, the projections of the 
eigenstates of H' n onto Vo, describing the kink at site n. Considering a 
time increment St such that the kink moves at site n + 1 while the coherent 
state representative point z changes to z + 5z = z + z5t, z §^ = 1, and using 
the defining properties of Vo, + <^)) turns out to be a superposition of 
{^ n+1 ' 1 }, describing the kink at site n + 1. Hence, for appropriate initial 
conditions, the kink moves spontaneously from site nton + 1, for all n. This 
proves that the system supports a kink-like excitation propagating along the 
chain at constant energy, once more as promised. 

In conclusion — We have proposed a quantum model for a 1-d chain of elec- 
trons and phonons, providing a self-consistent description of Peierls distor- 
tion and showing that the system is able to sustain soliton excitations. Key 
features of our approach are the non-linear coupling beteween electrons and 
phonons, which generalises the SSH model; the description of the phonon 
degrees of freedom in coherent-states, which permits a dynamical picture of 
the phonons; and the subsequent restoring of the original quantum dynam- 
ical symmetry. Further insight may be gained in the future via a thorough 
analysis of the ground state and of the soliton dynamics. This may open 
new perspectives for the application of solitons in one-dimensional quan- 
tum systems. For instance, solitons may be used as means of transferring 
quantum information. Work is in progress along these lines. 
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